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1.  Statement  of  the  Problem 

A  "hot-film  anemometer"  probe  Is  typically  sketched  on  Figure  1,  It 
consists  of  a  "probe  body"  made  of  an  Insulator,  with  a  very  thin  metallic  film 
deposited  on  Its  blunted  tip.  Two  conducting  wires  ("leadwires") ,  attached  to 
the  two  ends  of  the  film,  enable  the  latter  to  be  heated  electrically  by  a 
constant  external  d.c.  current  source.  If  the  probe  is  immersed  in  a  moving 
fluid,  then  the  cooling  effect  of  the  flow  lowers  the  film  temperature  and  hence 
its  resistance  and  voltage  drop.  A  detector  connected  in  the  circuit  can  then 
be  used  to  associate  steady-state  film  voltages  to  the  steady  state  flow 
properties,  or  the  a.c.  voltage  across  the  probe  to  fluctuations  in  these 
properties. 

The  essential  questions  are: 

•  (a)  How  Is  the  steady-state  film  voltage  (or  resistance)  quantitatively 
connected,  at  any  moment,  to  the  steady-state  flow  properties? 

(b)  How  is  the  film  a.c.  voltage  magnitude  quantitatively  connected  to  the 
magnitude  to  the  flow  fluctuations? 

(c)  How  Is  the  film  frequency  response  connected  to  the  frequency  of  the 
flow  fluctuations? 

These  questions  will  be  answered  in  this  paper  for  the  specific  case  of  the 
film  lying  on  the  stagnation  region  (point)  of  the  probe,  and  only  if  the 
electric  current  is  constant.  For  analyses  of  the  "constant  temperature" 
operation  (where  the  film  temperature  in  a  fluctuating  flow  Is  forced  to  stay 
constant  by  appropriate  variations  of  the  current)  and/or  different  probe 
geometries,  one  can  consult  the  literature.  For  example  Ling  (Reference  1)  and 
Ling  and  Hubbard  (Reference  2)  analyzed  the  general  case  of  the  present  problem, 
and  also  the  case  of  the  film  lying  on  a  30-deg.  wedge;  Lowel  and  Patton 
(Reference  3)  analyzed  the  case  of  the  film  painted  over  a  circular  cylinder; 
Bellhouse  and  Rassmussen  (Reference  4)  and  Blalvas  and  Zolkina  (Reference  5)  the 
case  of  the  constant  temperature  film  and  so  on.  These  various  solutions  are 
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generally  different  from  each  other  because  of  different  assumptions  regarding 
the  thermal  boundary  conditions  and/or  constraints  on  the  film  current  or 
temperature  (l.e.  constant-current  or  constant  temperature  operation).  The 
previous  work  closest  to  the  present  analysis  Is  that  of  Ling  and  Hubbard  in  so 
far  as  the  frequency  response  Is  concerned;  however  some  differences  exist  which 
are  fundamental  In  the  final  results. 

All  analyses  (see  Reference  6  for  some  early  bibliography)  Including  the 
present  one,  follow  the  same  general  path:  first,  the  power-balance  equation  for 
the  film  is  considered  jointly  with  the  conduction  equation  for  the  substrate, 
which  Is  frequently  assumed  one-dimensional.  The  film  equation  is  next  reduced 
to  an  ordinary,  first-order  non-linear  differential  equation,  which  then  is 
linearized  by  assuming  that  the  flow  fluctuations  are  small.  This  linearized 
equation  immediately  provides  (a)  a  steady-state  algebraic  equation  which  has  a 
simple  direct  solution,  (b)  an  equally  simple  ordinary  linear  first-order 
equation  for  the  film  temperature  fluctuations. 

Another  point  common  to  all  film  analyses  Is  their  resemblance  to  the  "hot 
wire  anemometer"  analysis.  The  latter  contains  a  different  type  of  conduction- 
loss  term  (Reference  7)  when  the  end  losses  are  considered,  and  no  conduction 
loss  term  at  all  for  the  "Ideal"  hot  wire;  and  this  Idealized  case  will  be 
frequently  mentioned  under  the  "hot  wire"  terminology  for  comparison  purposes  in 
this  report. 

2.  Basic  Film  Equations 

Neglecting  radiation,  and  referring  to  Figures  1  and  2  and  the  List  of 
Symbols,  the  film  power-balance  equation  Is: 

“  |T-T«>  ♦  Jsf+  *1) 

These  terms,  which  have  the  dimensions  of  power  flux  (energy  per  unit  area  and 
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time)  are,  from  left  to  right,  the  unsteady  film  temperature  change,  the 
convective  exchange  with  the  flow,  the  ohmic  heating  and  the  conductive  exchange 
between  film  and  substrate. 

The  convective  term  contains  the  Nusselt  number  defined  so  that  k0  is  the 
fluid  thermal  conductivity  based  on  the  flow  stagnation  temperature  T0,  related 
to  the  adiabatic  recovery  temperature  Te  via  the  recovery  factor  tj: 


T.  ■  i)T0 


The  characteristic  film  dimension  h  is  the  film  height  {Figures  1,2)  which  must 
then  also  be  used  to  form  the  film  Reynolds  number 

He,  -  (3) 


where  /i0  is  the  stagnation  viscosity.  For  steady-state  (ideal)  hot  wires  only 
the  first  two  terms  on  the  r.h.s.  of  (1)  are  present;  then,  since  in  that  case 
Af  -  AhM  -  *dl  and  h  becomes  the  wire  diameter  d,  the  hot-wire  counterpart  of  (1) 
will  be 


1ZR  -  *lk0N  (T-Te)  (4) 

In  the  case  of  films  Af  -  hw,  and  if  a  linear  resistance- temperature  relation  is 
assumed: 


R  -  Rr  (1  +  oT) 


(5) 
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Ra  «  Rr  U  +  «T.) 


(6) 


then  eq.  (1)  becomes 


e  g  -  -k.wN  (T-TJ  +  12Rr  (1+otT)  +  ksA,  ( 


(7) 


A  second  equation  entering  the  problem  deals  with  the  flow  of  heat  in  the 
substrate: 


=  B  (8) 

at  ax2 

where  the  film  is  assumed  infinitely  thin  so  that  T  -  Ts  (x-O). 

If  we  non-dimensional ize  by 

T  a  OtT,  Te  a  0T#>  T#  B  aT,,  t  B  X  a  x/JKf  (9) 

then  the  film  equation  (7)  Is 

,  3/2  _ 

SL.LRrO&f  (1+T)  -“WJi  (T-TJ+^f  (10) 

at  cB  eB  efl  dx 

and  the  substrate  equation  becomes 

i  -  &  (11) 

at  ax2 


with  boundary  conditions  to  be  discussed  below. 
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3.  Linearization  of  the  Film  Equation 

Sfnce  the  dependent  variable  T  of  the  problem  and  also  k0,  T,  and  N  are 
generally  time -dependent,  eq.  (10)  Is  non-linear.  Solutions  are  easiest  to 
obtain  when  all  these  quantities  exhibit  small  fluctuations,  denoted  by  a  prime, 
about  their  mean  (time-averaged)  values,  which  appear  unprimed,  as  follows: 


T  B 

0  +  0'  (t)  , 

,  d'  «  0 

(12) 

T.Bfl, 

+  (t)  . 

{13) 

T#  =  «B(x)  + 

(x,  t)  , 

(14) 

The  convective  term  coefficient  will  for  the  moment  be  written  as 

follows: 

+wAfX^l' 

■  f  +  f'(t) 

,  r  «  f 

(15) 

The  quantity 


B  b  (16) 

will  be  also  assumed  constant  (see  Section  10.4  for  exceptions).  If  equations 
( 12) - ( 16)  are  substituted  Into  (10)  and  products  of  two  or  more  small  quantities 
are  omitted,  we  obtain 

|^(d  +  *')  -  B  (1-rf+tf'J  -  f(tf-0e  +  *'-#',) 

-  f '(«  «.)  ♦  +  <#’•>«  1  <17> 

4.  Steady-State  Solution 

When  -  0'8  =  f'  -  0,  eq.  (17)  becomes 
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(18) 


b  (1+#)  -  f(tf  -  o.)  +  -  <> 

while  eq.  (11)  becomes 


^  -  0  (19) 

dir 

Eq.  (19)  can  supply  the  last  term  In  eq.  (IB)  If  two  boundary  conditions  are 
available  for  Its  solution.  One  of  these  conditions  is 


»a  (x  -  0)  -  «  (20) 

The  other  condition  is  specific  to  the  probe  geometry.  Often  the  substrate 
is  assumed  to  have  a  finite  thickness  with  its  far  face  at  Te,  but  this  condition 
is  sometimes  abused  by  misapplication  to  unsuitable  geometries  (Reference  4 
applies  this  condition  to  a  geometry  in  which  the  two  faces  of  the  substrate 
coincide).  In  the  present  instance  this  condition  Is  suitable,  as  Is  clear  from 
Figure  1.  The  probe  body  is  seen  to  have  a  large  length- to -diameter  ratio; 
furthermore,  the  wetted  film  surface  area  is  very  much  smaller  than  the  wetted 
area  of  the  probe  body.  Therefore,  if  the  film  is  electrically  heated,  the  probe 
body  temperature  will  rise  only  in  the  vicinity  of  the  film,  out  to  some  distance 
L  from  the  latter.  Beyond  that  the  probe  temperature  will  equal  its  adiabatic 
recovery  temperature  due  to  the  flow  and  It  Is  very  reasonable  to  assume  that  the 
film  and  probe  recovery  temperatures  are  equal  to  each  other  and  to  Te.  The 
calculation  of  L  will  be  discussed  In  the  next  Section;  for  the  moment  we  will 
assume  on  the  strength  of  the  above  arguments  that  the  second  boundary  condition 
Is 

T,  (x  ■  L)  *•  T#,  08  (x  -  L/yR,)  -  (21) 
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The  solution  of  (19)  subject  to  conditions  (20)  and  (21)  Is: 


f  (d#  -  d) 


(22) 


which  gives 


(“?)«« 

flx 


(23) 


Substituting  (23)  Into  (IB)  we  obtain  the  film  non-dimensional  temperature: 


5. 


*  +  - B - 

#  f  +  k#AVLcB  (24) 

fl  »  - n - 

1 - 0. - 

f  +  kBA*f/Lcfl 

Discussion  of  the  Steadv-State  Solution 

When  k,  -  0,  eqs.  (5),  (6)  can  be  used  to  simplify  eq.  (24)  Into: 


R - —  (25) 

i . 

iv*r 

This  is  the  steady-state  hot-wire  equation  which,  If  k0  Is  changed  to  ke  and  w 
absorbed  Into  N,  Is  Identical  to  eq.  (14)  of  Reference  8. 

With  the  substrate  present  (k,  *  0)  It  can  be  shown  that  (24)  reduces  to 


1  . 

k0wfl 

where  the  "compound  Nusselt  number" 


(26) 
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N  3  N  +  ^  A»  -  N  +  *«  (27) 
k0  wL  k0  L 

To  estimate  the  second  terra  on  the  right  of  (27),  which  is  the  steady-state 
"conduction  loss"  contribution  to  N,  we  need  to  estimate  the  "penetration  depth* 
L  first  discussed  in  the  previous  Section.  Anders  (Reference  9)  estimates  L  from 
a  simple  one-dimensional  steady  state  calculation  modeling  the  probe  as  a 
circular  cylinder  parallel  to  and  wetted  by  the  flow,  with  its  upstream  face  kept 
at  the  film  temperature.  He  finds,  setting  k,  -  k0,  that 


L 


[k*rpLe 

20“ 


1/2 

.] 


o  p 


(28) 


If  the  probe  characteristic  length  Lc  Is  identified  with  the  cylinder  (probe 
body)  radius  rp,  then 


N  =  N  +  ^  [2Np  (J*)]VZ  (29) 

For  typical  materials  and  conditions  the  convection-loss  term  on  the  right  of 
(29)  is  of  order  1-10,  and  may  be  larger  than  N  itself. 

For  hot  wires  the  conduction  loss  formula  due  to  the  wire  supports  (see 
Reference  10  for  discussion  and  additional  References)  which  corresponds  to  (29) 

is 

N  *  N  +  4  d*e  ** 

l2  K0 

with  typical  values  of  order  0.1  for  the  second  term  on  the  right,  the  only 
similarity  with  the  second  term  of  (29)  being  the  appearance  of  the  same  type  of 
conduction-loss  parameter  (a  ratio  of  lengths  and  one  of  conductivities).  This 
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parameter  appears  In  the  square  root  In  eq.  (29),  and  In  the  first  power  in  (30), 
and  also  the  film  loss  term  contains  k,  while  that  for  the  wire  contains  The 
principal  difference,  however,  is  one  of  magnitude:  based  on  the  numbers  quoted 
the  conduction  loss  for  the  hot  wire  is  a  mere  correction  to  its  measured  value, 
while  for  the  film  the  correction  overwhelms  the  Nusselt  number  Itself. 

6.  Unsteady  Solution 

6.1.  The  Time -Dependent  Equation 

When  eq.  (IB)  is  subtracted  from  (17),  the  unsteady  equation  for  the  film 
temperature  perturbation  6'  is  obtained: 

in 

*  «'  (f-B)  -  -  f '  (<-«,)  (31) 

and  the  substrate  equation  becomes 


«*'.  .  **•'.  (32) 
ST-  3x* 

The  r.h.s.  of  eq.  (31)  is  the  forcing  function  driving  the  linear  system 
on  the  left.  In  principle  f'  and  0' t  are  caused  by  flow  fluctuations  alone  (cf. 
eqs.  (13)  and  (IS)).  This  will  be  shown  in  a  later  Section,  but  an  important 
exception  requires  discussion  at  this  early  stage.  Since  f'  Is  driven  by  N',  N' 
by  r'  and  r'  by  #'e  and  a  term  containing  9'  will  appear  on  the  right  of  eq. 
(31).  The  remainder  on  the  r.h.s.  will  depend  only  on  the  flow  fluctuations, 
assumed  stationary  and  thus  expressible  as  a  Fourier  series.  For  a  moment  we  are 
concerned  with  one  of  the  Fourier  components,  of  frequency  «.  The  amplitude  of 
this  component  must  have  the  form  E/t  where  r  Is  the  time  constant  of  the  system 
(this  procedure  is  standard  for  "first  order"  systems,  which  In  this  case  means 
the  no-substrate  case).  The  quantity  E  will  Include  the  amplitude  of  the  flow 
fluctuations,  and  its  calculation  will  not  be  needed  until  Section  7. 
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Collecting  the  ideas  expressed  In  the  above  paragraph  ,the  flow  forcing  terra 
(as  can  be  seen  considering  eq.  (75J  to  be  presented  later)  is 


«'e 


f '(#-«,) 


I 


si  nut  -  fNr#' 


(33) 


so  that  eq.  (31)  Is 


.3/2 

#'[f(l+Hr)-B]  -  sinwt  (34) 


The  boundary  condition  for  this  equation  is  to  set 

0'  -  A  $in(5t  -  f)  (35) 

where  A,  $  are  to  be  determined. 

6.2  Unsteady  Solution  in  the  Absence  of  the  Substrate 

We  first  consider  the  absence  of  the  substrate  (k,  *»  0),  In  which  case  eq, 
(34)  becomes 


$-+  r  [f(l+Nr)-B]  -  |s1n(K)  (36) 

This  is  the  familiar  linear  ordinary  differential  equation  of  the  first  order 
with  constant  coefficients,  which  describes  the  so-called  "first  order"  systems. 
The  inverse  of  the  second  term  coefficient  Is  the  (non-dimensional)  system  time 
constant: 

if  ■  -nrabr  (3?) 

The  solution  of  (36),  given  in  all  standard  texts,  Is 
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- -  sin  (it  -  i) 

[l+(«r)2]1/z 

where  £  b  tan*1  ur.  Note  that  «  *  fi^f>  so  that  it  =  wt. 


(38) 


6.3  Discussion  of  the  Unsteady  Solution  Without  Substrate 

If  we  recall  the  definitions  of  f  and  B  from  equations  (15)  and  (16),  the 
dimensionless  time  constant  given  in  (37)  reduces  to 


Afk#Nw(  |j.  +  Nr) 

and  Its  dimensional  counterpart  to 


(39) 


wk0N  (g«  +  Nr) 


(40) 


The  time  constant  t  given  in  (39)  Is  the  "inherent"  film  time  constant,  as  if  the 
film  was  "suspended  in  space”  without  a  substrate.  For  this  reason  it  is  also 
the  "hot  wire"  time  constant,  familiar  to  most  readers  in  the  case  Nr  -  0,  for 
which  we  obtain  the  classical  formulas  (e.g.  Reference  7,  p.  241): 


t  -  — -  €  (lt£i  ,  r  (r=0)  =  (41) 

wk0NRe  wk0N  wk0N 

and 

— - —  =  1  +  r  (42) 

7  (1-0) 

Note  that  in  practical  cases  these  formulas  differ  only  slightly  from  the  exact 
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expression  (40),  since  Nr  «  1  usually. 

The  bracketed  factor  in  the  denominator  of  (38)  Is  the  familiar  attenuation 
due  to  the  inherent  time  constant  of  the  hot  wire  (or  the  no-substrate  film). 
It  Is  seen  that  at  large  u  the  signal  E  is  attenuated  by  8  db/octave  {a  factor 
10  In  signal  drop  for  a  factor  10  Increase  in  frequency);  see  Figures  3  and  4. 

6.4  Solution  of  the  Unsteady  Problem  with  Substrate 

To  solve  the  complete  equation  (34)  when  k,  t  0,  eq.  (32)  must  first  be 
solved  for  #'s.  The  boundary  conditions  for  (32)  are: 


t,  t  s  0  :  »  0 

(43) 

t,  t  >  0  :  fl'a  (x-0)  -  6' 

(44) 

In  view  of  eq.  (35),  this  last  condition  is  rewritten  as: 

d'a  (x-0)  -  A  sln(ut-d)  (45) 

The  solution  of  eq.  (32)  with  conditions  (43)  and  (45)  Is  a  well-known  one  (e.g. 
Reference  11,  p.  93  or  Reference  12,  p.  112)  and  is: 

0'a  (x,  t)  -  Ae  2  sin(wT-d-xV^)  (46) 

and  ( jf- )**  -  -A/^tsIn  (ut  -  d)  +  cos(«t  -  *)]  (47) 

With  the  aid  of  (35),  (37)  and  (47),  eq.  (34)  becomes 
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Awcos(wt  -  4)  '*  f  sin  (®t  -  0) 

3/2  _  p,  - 

+  A/^  [sin  (tit  -  tf)  +  coslttt  -  d)]  “  f  si  nut  (48) 

By  zeroing  the  coefficients  of  sinwt  and  cosut  we  solve  for  A  and  0  for  Insertion 
into  eq.  (35);  the  result  for  the  non-dimensional  film  perturbation  d'  and  phase 
angle  ^  is: 

fl'(t)  -  o  T'(t)  - - E  slntikA) -  (49) 

t(«f+qvSf)z  +  (i+q^f)z],/2 

i  -  tan'1  t  fl. Jjf  ■  (50) 

i  +  qyuf 

i 

where  f  is  as  already  given  in  (39)  and  Q,  which  will  be  called  the  "conduction 
loss  factor"  or  simply  "loss  factor"  Is 


Q  . 


fB 


(51) 


An  alternative  form  of  Q,  useful  for  practical  computations,  is  given  in  Section 

10.2. 


6.5  Mscussion  of  the  Unsteady  Solution  with  Substrate 
Since  ttr  *  t)7 1  the  "attenuation  factor" 


£ 

E 


(52) 
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and  the  phase  angle  ^  are  plotted  on  Figures  3  and  4  for  a  range  of  Q.  When  q 
-  0  there  is  no  substrate,  and  the  attenuation  Is  of  the  hot  wire  type  (l.e.  eq. 
(38)).  Notice  that  as  soon  as  substrate  Is  added  to  a  specific  film,  A/E  always 
drops  l.e.  the  attenuation  Increases.  This  Is  a  vital  point,  since  the  question 
Is  often  asked  whether  "hot  films  are  better  than  hot  wires".  This  question  is 
wrong  because  it  Is  too  general.  If  a  wire  Is  compared  with  a  film  of  the  same 
Inherent  time  constant  f,  then  the  film  Is  always  worse  (l.e.  It  attenuates  more) 
than  the  hot  wire.  This  Is  true  for  a  certain  frequency  range  only;  at  very  high 
frequencies  film  and  wire  are  "Identical",  l.e.  their  attenuation  is  of  the  hot 
wire  type  (6  db/octave). 

The  3  db/octave  film  attenuation  characteristic  was  originally  mentioned 
by  Ling  (Reference  1)  and  Ling  and  Hubbard  (Reference  2).  A  graph  published  In 
the  latter  Reference  compares  the  film  to  the  hot  wire  quite  favorably  but  also 
improperly,  because  there  is  no  accompanying  mention  of  the  Inherent  time 
constant  and  loss  factors. 

7.  Film  Probe  Response  to  Specific  Flow  FI uctuatlQliS 
7.1  Introductory 

In  the  foregoing,  the  various  types  of  flow  fluctuations  (l.e.  velocity, 
density  etc.)  have  been  Indiscriminately  represented  by  the  factor  E  in  equation 
(49),  so  as  to  first  focus  attention  on  the  frequency  response  Issue.  The  next 
objective  is  to  express  E  In  terms  of  these  fluctuations,  l.e.  to  analyze  the 
"forcing  function”  appearing  on  the  r.h.s.  of  eq.  (31).  Specifically  we  want  to 
put  the  forcing  function  Into  the  form 


where  v,' 


-  f'  («-«.)  -  2  e, 

Is  the  time-dependent  small  fluctuation  In  some  flow  property 


(53) 

(e.g. 
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velocity  u'  or  density  p'),  v,  the  corresponding  mean  value  and  e,  the 
corresponding  sensitivity  coefficient. 

At  the  outset,  it  must  be  said  that  since  the  forcing  function  is  Identical 
to  that  of  a  hot  wire,  the  answer  must  be  identical  to  that  supplied  by  Horkovin 
(Reference  13).  Here  we  will  therefore  supply  results  which  are  the  same  as  - 
Morkovln's  with  the  following  differences:  (a)  No  circuit  Impedance  losses  are 
considered,  (b)  the  controllable  variable  is  the  overheat  r  rather  than 
Morkovin's  A'H,  (c)  the  resistance-temperature  relation  is  herein  assumed 
1 1near. 

7.2  Fundamental  Relations  for  Small  Fluctuations 
Since 


-  f'(<M>e).  -  f  (*-*.)  [  -  7-] .  <54> 

f 

we  will  first  consider  the  dependence  of  the  bracket  on  the  fluctuations 

V  K<  \  (total  temperature)  (55} 

u'  «  u  (velocity)  (56) 

p'  «  p  (density)  (57) 


where  the  r.h.s.  of  these  inequalities,  without  primes,  represent  the  mean 
values.  Assume,  with  Morkovin,  that  the  Nusselt  no.  N  depends  on  Mach  number  M, 
overheat  r  =  (R  -  Re)/R#  and  Reynolds  number  based  on  stagnation  conditions: 

Re  =  ^  (56) 


so  that 
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N  >  N  (H,  Re,  r) 


(59) 


Small  fractional  fluctuations  (denoted  by  either  a  differential  or  a  prime)  then 
are: 

IP  ■  IF-  -  J  (ft,  «•  ♦  SB  dH  .  fjj  dr)  . 

where,  for  brevity,  we  symbolize  for  any  two  quantities  P  and  Q: 

p  H  Q  SS. 

K0  H  P  3Q 

Now  assume  the  viscosity-temperature  relation 

M  -  Tn,  Mo  -  V 

which  gives 

r-  S“  *  ■  IF+  f-"lr 


ttl+  N  El 

H  +  Nr  r 

(60) 

P  and  Q: 

(61) 

(62) 

/ 

(63) 

For  dM/M  we  utilize  Morkovin's  expression  (eq.  (25)  of  Ref.  13): 

“B-B1** <iH  lr>;  * 


(64) 


Since,  upon  assuming  a  reference  temperature  Tr  *>  0, 


R  =  Rr  (1  +  oT)  ,  R.  -  Rr  (1  +  ®Te) 


(65) 


and 


.B.-.R,  .  aOJoi 

R,  1  +  arta 


(66) 
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we  also  get 


dr 

r" 


r' 

r~ 


V  T,(1«T)  T'e 
T^r;  -  (T-Te)  ll+al,J  TT~ 


and  since  T,  =  jjT,: 


V 


nl 

*1 


The  recovery  factor  dependence  for  air  is  known  to  be 


n 


*(H,  Re),  *1-5 
t]  t] 


°  % 


Re' 

“Re 


Combining  equations  (63),  (64),  (66),  (67)  and  (68), 


rf 

r~ 


Te(3+aT) 
IT-1,) (1+alJ 


♦  %  *(  -- 


1  V 


) 


Equations  (63),  (64)  and  (70)  can  now  be  substituted  Into  (60): 


N'  HdN 

r  r 


'Re 


(u'+^  - 
U  fl 


)  +  Nk  6  (5  - 
u 


1 

"2 


.  NrT,(l+oT) 
(T-T.)(1«T.) 


+  %  «(  _U'* 


1 

"2 


] 


In  addition,  we  need  the  fluctuation  In  the  thermal  conductivity,  which 
consider  to  be 


(67) 


(68) 


(69) 


(70) 


(71) 
we  can 


(72) 
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so  that 


*■'  -  m  V  {73) 

which  concludes  the  collection  of  fundamental  relations  needed. 

7.3  Development  of  the  Forcing  Function 

He  now  develop  the  bracket  of  eq.  (54)  with  the  aid  of  the  definitions  of 
f»  (eqs.  (12),  (13)  and  (15))  and  of  relations  (63),  (64)  etc.: 

V  f'  f'  -  Te  Te'  (  K'  +  N') 

n;  "  r  r  t;  v  it 

(74) 

Expanding  the  fluctuations  of  k0,  Te  and  N  into  their  constituents, 


Cl  ♦  %  * 

V  u 


1  T0 
T0 


)  +  %.( l'+ 

U 


-m  V 


)  +  N„  «(“' 
u 


1  V)  +  NrT'  .NrT,(l+aT) 
U  ^  (T-Ta)  (  1-hjT,) 


{]°-  +  *ln&  I  It  )  +  *«.(  —  *  fl  -  n  I?  ))] 

*  u  7  TT  u  p  in 


(It  W'tl**T>  )  (!-«%. 
l+ar. 


®  +  +  ^  Nm  -  m) 
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*“1  tJu  (1*  W'(1*°T1  )  l«*  +  <k.\  '  -  «N,] 


*  *1  J.  (1+MWH  V' 


T+aT 


(75) 

The  last  term  In  this  equation,  containing  T',  will  be  omitted  from  further 
considerations  since  it  has  already  been  carried  to  the  l.h.s.  of  eq.  (34).  (See 
discussion  in  Section  6.1).  The  complete  forcing  function  as  given  by  (54)  can 
now  be  found  by  noting  that  the  factor 


f  ffl-0  1  -  MoN"  o(T-T  1  -  “  r  R« 

'  TB  •  Tr 

To  get  the  r.h.s.  of  eq.  (34)  we  now  multiply  this  factor  with  eq.  (75): 

W-V  {JC  -  f'  )  -  A<kq.Nw  R*  (Til [G  -  \ 

0-0,  f  fB  Rr  T0  c 

+r(nNBe  +  _£  NH  -">)]+  Ji'tG  (H,  +  W  *  r(N„e  +  «*,)] 


(76) 


+  L  -  ■*„.])  <77> 

/> 

where  the  T'-containlng  term  of  eq^  (75)  is  missing  by  design  .as  noted,  where 

G  a  aT«  (1  +»r(l+ttT)  )  ,  ftT»  fl+N.(r+l)]  (78) 

1+aT,  l+ale  l+aTe“ 


and  where  t}H ,  N„e  etc.  have  been  defined  via  eq.  (61)  and  $  via  eq.  (64). 


Eq.  (77)  can  be  furthermore  re-written  If  we  utilize  the  definition  of  the 
non-dimensional  film  time  constant  T  from  eq.  (39)  as  follows: 
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f(*‘*.)  ( 


Afk0Nw  R. 


cB  Rr 


"r*  fr 

»  *fr 


Afk0Nw  Re 


{ . )  - 

*  [Nr  *  +  lr1  {....) 
K  ^ 


(79) 


and  eq.  (77)  can  then  be  rewritten  with  the  factor  outside  the  curly  bracket 
replaced  by  (1/f) <R/Rr)  and  with  G  and  r  replaced  by 


G  a 

(80) 

T+al7 

r  ■ 

r 

(81) 

r-  r 

(r+1)  Nr+1 


8.  Final  Solution  for  the  Film  Temperature  Perturbation 

The  fluctuations  To'/To,  u'/u,and  p'/p  in  eq.  (77)  have  been  assumed  to  be 
functions  of  the  time.  If  each  Is  representable  by  a  Fourier  series  In  sln(ut) 
and  if  we  concentrate  on  the  Fourier  component  at  a  specific  frequency  a,  we  can 
write 

V  b  t  slnat,  ^  -  u  slnut,  -  p  slnut  (82) 

T;  up 

where  (  )  stands  for  the  non-dimensional  amplitude  (e.g.  T0  Is  the  small 
fluctuation  amplitude  In  Ta,  divided  by  the  mean  value  of  T0).  In  view  of  eq. 
(82),  we  can  then  re-write  eq.  (77)  as  follows: 

f  «-*.)  (JL*.  ■  »  >  -  i  "  (T.  [S(l-ni|„  -  \) 
f  T  Rr  Z 

♦  ?  (nil,,  ♦  |  n,  -  ■)]  t  5  [6  (i^  t  m,.)  -  f  (n,„  +  <«,)]  (B3) 

*  foil,.  ■  ?"„]) 

But  we  note  that,  since  ur  -  ur,  this  forcing  function  has  exactly  the  same  form 
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as  that  on  the  r.h.s.  of  eq.  (34),  with 


E 


(84) 


Therefore,  since  the  solution  for  eq.  (34)  Is  given  by  eq.  (49),  by  analogy  we 
see  that  in  the  present  Instance,  the  solution  for  fl'  Is 


aT'  ■  d'(t) 


R  _ slnf«t-*l  - 

Rr  [  (wt  +  QVSZ  +  (1  + 


(Tp  [G  d-n»7,e 


-  *  »„)  +  r  (nN„.  +  *  Nm  -  m)3  +  G  [G  ( &vH  +  »Re) 

-  r  (HR#  +  6N„)]  +  fl  tGq#e  -  rNRe]}  (85) 

and  with  0  and  Q  remaining  as  given  In  eqs.  (49)  and  (50)  and  r  or  f  as  given  in 
(39)  and  (40).  (Note  again  that  uf  -  «t). 


9.  Film  Voltage  Fluctuations 

The  fluctuation  e'  in  the  film  voltage  can  be  obtained  directly  from  eq. 
(B5),  since  It  is  related  to  the  mean  film  voltage  e  as  follows.  Since  for 
constant  current 

e  +  e'  =  1(R  +  R'}  «*  1Rr  [1  +  o  (T  +  T')]  (86) 
then  e'  ■  1aRrT'  ,  e  ■  iR  (87) 
So  that  by  multiplying  both  sides  of  (B5)  by  iRr  we  get 


jRsio(rt-tf)  eglnUrt-#)  . 

6f  =  ioRrT'  »  (•>>•)  *  {>•■■} 

(  }m  [  11/2 


10.  Practical  Aspects  of  the  use  of  the  Film  Anemometer 
10.1  Mean-Flow  Measurement  Routine 


(88) 
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The  universal  mean  heat-transfer  characteristics  of  films  (i.e.  the  N  (Re, 
M,  r)  and  q  (Re,  H)  dependences)  are  of  fundamental  Interest  only  If  they  are 
composed  from,  and  afterward  applied  to,  films  of  accurately  known  geometry. 
Such  accuracy  is  impossible  for  home-made  film  probes*,  and  universal 
characteristic  curves  replacing  calibration  curves  will  therefore  not  be 
available.  However,  a  specific  film  probe  can  be  calibrated  and  then  used  for 
a  measurement  without  the  need  of  knowing  Its  exact  geometry.  In  all  cases, a 
measurement  of  the  Nusselt  number  Is  needed.  The  procedure  for  this  measurement 
Is  generally  Indicated  in  Reference  8. 


According  to  eq.  (27), 


N  -  N  +  k*  h 


or 


wN 


wN  + 


k8  hw 

^  "C" 


(89) 


(90) 


According  to  Reference  8,  what  Is  usually  measured  (and  plotted  vs  Re'  or  Reh) 
Is  the  wN  on  the  l.h.s.  of  eq.  (90).  (In  that  Reference  the  symbol  N  stood  for 
wN  as  defined  here).  If  w  Is  known,  N  can  next  be  obtained  from  eq.  (89) 
provided  that  (kg/k0)  (h/L)  can  be  computed,  which  is  Improbable  because  of 
uncertainties  In  L.  However,  (kg/kj  (h/L)  does  not  depend  on  Re',  and  can  be 
obtained  from  the  Intercept  at  Re'  ■  0  of  a  plot  of  N  vs  Re'.  Thus  a  plot  of  N 
vs  Re'  is  obtainable. 


*  By  contrast  hot  wire  probes,  even  of  the  home-made  variety,  enjoy  the 
precision  established  by  industrial  standards  for  the  wire  material,  e.g.  Its 
diameter  and  resistivity. 
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.  If  w  is  unknown  eq.  (90)  is  operative,  the  Intercept  just  mentioned  will 
signify  (k,/k0)  (hw/L)  rather  than  (k^kj  (h/L),  and  the  final  result  will  be  a 
value  of  {wN)  for  each  Re'  or  Reh. 

When  the  final  objective  is  to  deduce  Re'  in  an  unknown  flow  by  measuring 
the  Nusselt  number  of  a  specific  film,  it  is  adequate  to  use  a  film  calibrated 
as  (wN)  vs  Re',  provided  that  the  total  temperature  for  calibration  and 
measurement  Is  the  same  (this  will  ensure  that  (k,/kft)  (wh/L)  Is  the  same).  If 
the  two  total  temperatures  differ  by  a  known  amount,  changes  in  the  conduction 
loss  can  still  be  estimated  based  on  changes  In  k0. 

For  fluctuation  measurements,  it  has  been  seen  In  Section  7.2  that  the 
Nusselt  no.  enters  the  sensitivities  via  its  logarithmic  derivatives,  e.g. 

N  Re  dN  Re'  dwN  dwN  Re'  (91) 

■  —  —  “  -  -  °  -  — = - c — rx 

N  flRe  wN  5Re'  aRe'  wN  -  ™ 

which  can  be  computed  from  a  knowledge  of  the  wN(Re')  curve  and  the  determination 

of  the  conduction  loss  from  that  curve. 

10.2  The  Unsteady  Loss  Factor  0  Due  to  the  Substrate 

In  eq.  (49),  It  was  shown  that  the  film  probe  attenuates  unsteady  signals 
in  a  way  dependent  on  two  parameters:  the  inherent  film  time  constant  7  and  a 
loss  factor  Q.  For  finite  Q  the  film  response  is  quite  complicated  (see  Fig.  3) 
and  one  cannot  use  electronic  compensation  circuits  employed  for  hot  wires. 
Therefore  some  knowledge  of  the  magnitude  of  Q  is  needed  for  typical  films. 

The  loss  factor  Q  for  zero  overheat  (unheated  film)  is  given  by  eq.  (51) 

which  can  be  re-written  in  the  following  way« 
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3/2  rr 

Q  a  C“  1 1 

r  r  a 


{^) 

/»fct 


(92) 


To  make  a  typical  calculation,  assume  that  the  substrate  (subscript  s)  Is  quartz 
glass  with  the  following  properties  (which  are  taken  from  Reference  14  for  0°  C, 
with  the  page  numbers  noted). 

-  Quartz  thermal  conductivity  (p.  1867):  k,  -  0.023  cra  secldeg  C 

(93) 

-  Quartz  density  (p.  1689):  p,  -  2.65 

(94) 

-  Quartz  heat  capacity  (p.  1790):  c8  -  0.235 

(95) 

The  ke  is  the  thermal  conductivity  of  stagnant  air  at  0°  C  on  the  surface 
of  the  unheated  film: 


k*  "  x  1°  cm. sec  deg  C 

(96) 

The  film  material  Is  taken  to  be  platinum,  and  although  ultra-thin  films 
usually  have  properties  differing  from  those  of  bulk  materials  (Reference  6), 
bulk-material  platinum  properties  will  be  assumed,  again  at  0  deg.  C.: 

-  Platinum  density  (p.  461):  pt  -  21.45  gm/cm5  (97) 

-  Platinum  heat  capacity  (p.  1773):  cf  =  0.03162  gf/jjL  C 

(98) 

The  film  dimensions  will  be  assumed  to  be 

h  -  0.04  cm  (99) 

w  =  0.1  cm 
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As  Ling  points  out  (Reference  1),  the  film  thickness  d  can  be  computed  from 
Its  resistance,  using  a  resistivity  of  platinum  (at  0°  C)  of 

rK  -  10.96  mlcroohm.cm 

We  will  assume  as  typical  a  film  resistance  R  ■  20  ohms  at  O^C;  then 

d  ■  filf  - 1,3  x  10"6  cm  (100) 

Finally,  we  will  adopt  values  of  the  Nusselt  number  N  taken  from  measurements  of 
the  compound  Nusselt  number  N  already  made  at  MSU  (Reference  8);  for  example, 

N  -  30  (101) 

If  the  values  (93)  -  (101)  are  substituted  into  (92),  one  obtains 

Q  -  436  (102) 

This  loss  factor  means  that  the  substrate  completely  dominates  the  film 

performance;  see  Figure  3. 

10.3  The  Inherent  Time  Constant  of  the  Film 

To  complete  the  computation  of  attehuation,  the  value  of  r  (the  inherent 
time  constant)  is  needed.  This  can  be  obtained  from  eq.  (40).  Here  we  compute 
the  unheated  inherent  r  (for  r  =  0),  using  the  numerical  values  given  by  (93)  - 
(101)  with  added  assumption  that  Nr  =  0: 

r  -  e  -  *fcfhd  -  21  x  10'6  sec  (103) 

W  i gr- 

Flgure  5  shows  that  with  the  substrate  in  place  the  attenuation  is  much  larger 
than  without  it.  This  film  attenuates  more  than  even  a  large-time-constant  hot¬ 
wire,  e.g.  one  of  200  microseconds. 

10.4  Measurement  of  the  Response 
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Because  of  the  sensitivity  of  the  frequency  response  to  the  local  flow 
conditions  {via  the  l^N  or  kgN  parameters),  any  response  measurements  needed  to 
verify  eq.  (52)  must  be  made  "in  situ".  The  most  convenient  such  measurement  Is 
to  record  the  attenuation  caused  by  the  film  to  an  Impressed  time-dependent  (e.g. 
periodic)  heating. current.  The  issue  arising  here  Is  whether  the  attenuation  Is 
the  same  for  such  a  variable  current  as  It  would  be  for  variable  flow  properties. 
Assurance  to  this  effect  Is  mentioned  by  Kovasznay  for  hot-wires  (Reference  7) 
and  by  Ling  for  films  (Reference  1).  Proof  of  these  assertions  will  now  be 
given. 


10.4.1  Checking  the  Response  When  the  Substrate  is  Absent 

First  we  Investigate  the  response  of  the  film,  without  substrate,  Immersed 
in  a  quiescent  flow: 

T0'  ■  u'  ■  p'  ■  0 

So  that  the  forcing  function  on  the  r.h.s.  of  eq.  (33)  is  simply* 

-  f'(#-de)  -  -fNr*'  (104) 


when  the  heating  current  is 

1  +  i'(t),  1  steady,  1'  «  1  (105) 

Assume,  furthermore,  the  harmonic  input 

i'/1  *■  I  slnwt,  I  «  1  (a  pure  number)  (106) 

Inserting  (104)  and  (105)  into  the  general  equation  (17)  for  small  perturbations, 
and  considering  the  substrate  absent,  we  get 


*For  example,  consider  the  r.h.s.  of  eq.  (75)  when  To'  *>  u'  -  p'  -  0. 
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d«'  .  aR^f  C1z  +  211')  (1  +  0  +  0')  -  f(0  +  0'  -  0e)  -  fNr0'  {107} 
3T  ~7B 

which  in  the  steady  state  becomes 

B  (1+0)  -  f(0-0e)  -  0  (108) 

giving  eq.  (24)  as  the  final  sveady-state  result.  Next,  the  unsteady  equation 
is  obtained  by  combining  (107)  and  (108): 


{Jf-  +  0'[f(l  +  Nr)-B]  -  2B  (1+0)  (109) 

The  r.h.s.  of  (109)  is  a  forcing  function  replacing  the  forcing  function  on  the 
r.h.s.  of  (31),  when  k,  -  0.  Equations  (31)  and  (109)  are  Identical  otherwise, 
and  therefore  their  time  constants  are  the  same,  proving  Kovasznay's  statement. 
Therefore  the  film  time  constant  found  by  Imposing  a  current  fluctuation  Is. the 
same  as  that  prevailing  under  flow  fluctuations  in  the  absence  of  substrate. 

It  is  useful  to  proceed  with  the  solution  of  (109)  which,  according  to  eq. 

(38): 


sin(wt  -  0) 

-  orT'  *  0'{t)  -  [2B(l+0)  r  I]  <110> 

with  f  from  (39).  With  the  aid  of  eq.  (40)  the  bracketed  factor  in  (110)  is 

2B(1+0)?I  =  21iZ“rZ  (111) 

Now  note  that  since  the  total  voltage  drop  across  the  film  is  e  +  e'(t)  =  1R  + 
e'(t),  we  have 
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e'  =  d(  1 R)  -  Ri'  +  1oRfT'  (112) 

and  therefore,  utilizing  eq.  (110), 

-  Islnwt  +  2Ii2qR2  (113) 

e  R  k0NR.  (1  +  tV)1/2 

«  I  (slnwt  +  2r  s1ntft>t"^  )  (114) 

e  (l+rV),/z 

In  many  practical  circuits  the  Islnwt  component  of  e'/e  Is  electronically 
cancelled  ("resistive  component  cancellation"),  giving  a  "net"  film  voltage 
signal 


e'„  -  2rl  sin(wt'^) 

7~  (lirV)w  (115) 

This  equation  can  be  used  for  one  of  the  many  possible  ways  of  measuring  r.  For 
example,  the  mean  square  (m.s.)  of  e'n  at  a  certain  frequency  a  can  be  compared 
with  Its  equivalent  at  o  -  0;  then 


j  »  1  m.s.  tfil  -  ]]'^z 

“  e'  (u) 

1,0.4. 2  Checking  the  Response"#^  the  Substrate  Present 


(116) 


It  remains  to  consider  the  effect  of  a  current  fluctuation  In  the  presence 
of  the  substrate.  The  development  follows  closely  the  lines  of  the  previous 
section.  We  begin  again  by  considering  the  effect  of  equations  (104)  and  (105) 
on  eq.  (17): 


dfl'  c  aRrAfi2  +  Bd,  +  Zii'aRrAf  (i+fl)  -  f(#-#J 

ar  eB  IF- 


3/2 

^F  3F>-° 


+  k*A'  [(Ms) 


♦  (fiij  1  •  •'«, 

33T*=o 


(117) 
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In  the  steady  state  this  reduces  to  eq.  (18)  once  again,  with  the  solution  given 
by  (24).  When  (18)  and  (117)  are  combined,  the  unsteady  equation  is- 

S/2 

W'  +  0'  [f(l+M.)-8]  -  k.  (0tf'*)  -  2  1#  B  (1+*')  (118) 

3T  r  1  7B¥w  T 

This  equation  is  Identical  to  eq.  (31)  except  that  the  forcing  function  is  again 
now  driven  by  the  fluctuating  current.  Therefore  the  time  constants  of  the 
systems  described  by  eqs.  (31)  and  (118)  are  the  same,  proving  Ling's 
proposition. 


We  conclude  with  the  solution  of  (118).  If  we  take 

0'  *  A  sin  (wt-0)  (119) 


V 

1 


Is i nut 


(120) 


we  again  obtain  eq. 
Au  cos(wt-d)  +  zr 

T 


147),  giving  the  following  result:  ' 

3/2  _ 

sin(ut  -  $)  +  k*Af  A  [sin(ut-d)  +  cos(«t-0)] 

nr  z 


»  2IB  (1+0)  sinut 


(121) 


The  amplitude  A  and  phase  angle  d  are  found  once  more  by  zeroing  the  coefficients 
of  cosut  and  slnut,  and  by  utilizing  eqs.  (Ill)  and  (112)  we  find  that  in  this 
case  the  counterpart  of  eq.  (114)  is: 


f'.  Ksinut  +  2r _ sln<Mt+> _  )  (122) 

e  I(«r  +  qySr)*  +  (1  +  Q  Jut)1]'12 

The  net  r.m.s.  e'n  of  the  voltage  across  the  film  (i.e.  when  the  resistive 
component  Isinwt  is  cancelled)  is 


e 


IrJT 


(123) 
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The  objective,  which  Is  the  knowledge  of  the  bracket  In  the  denominator  of  (123), 
can  now  be  attained  by  measuring  all  other  properties  in  eq.  (123)  at  any  moment 
of  the  film  operation  In  a  given  flow. 

Finally,  It  may  appear  somewhat  odd  that  the  voltage  fluctuations  given  in 
eqs.  (115)  and  (122)  depend  on  the  overheat  r,  especially  since  this  would  make 
the  time-constant  measurement  difficult  at  low  r.  The  explanation  Is  that  the 
impressed  current  fluctuations  1'  are  assumed  to  be  a  small  fraction  of  the 
steady-state  1.  Therefore, as  1  (and  hence  r)  tends  to  zero  so  does  the  assumed 
i'  and  hence  also  the  voltage  fluctuations  e'n.  It  would  be  valuable  to  analyze 
the  case  where  1'  Is  small  but  nevertheless  much  larger  than  any  steady  current. 
This  would  represent  the  "a. c. -heated"  film. 


The  sole  reason  for  the  usefulness  of  the  hot-wire  at  high  speeds  is  the 
simplicity  of  Its  attenuation  characteristic,  given  by  eq.  (38): 


eour  _ | eouT 

efn  [1  +  uV]1'2  ur  -  -  eln  u 


(124) 


for  which  a  compensating  amplifier  with  a  gain  Increasing  as  u  can  be  built. 


To  equip  the  film  probe  with  a  similar  compensator,  one  would  hope  for 
similar  simplicity  In  its  attenuation.  Pessimism  In  that  direction  prevails  in 
the  earlier  hot-film  analyses,  which  replace  eq.  (124)  by  complicated  functions 
of  u,  due  mainly  to  boundary  conditions  different  from  those  used  here.  For 
example,  In  References  3,  4, and  5, one  finds  attenuation  functions  with  minima, 
maxima,  regions  of  insensitivity  to  u.etc.  We  shall  now  show  that  the  present 
results  simplify  this  function  considerably,  with  hopes  of  building  its 
electronic  Inverse  for  compensation  purposes. 
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The  attenuation  found  herein  is  given  by  eq.  (49): 

eouT  _  _ 1 _  (125) 

ein  E(«r  +  q^5r)z  +  (1  +  Q,/ut)2)1/z 

This  function  Is  not  necessarily  a  smooth  one  In  u  since  it  goes  as  the  Inverse 
square  root  of  frequency  up  to  some  w,  and  as  the  Inverse  frequency  thereafter. 
However,  if 


q  »  1  (126) 

Q  »  jm  (127) 

then  in  our  present  case  eq.  (125)  gives 

ecuT  _L_  (128) 

V  “ 

In  Section  10.2  we  quote  estimates  of  Q  -  436,  r  *  20  microseconds.  For  these 
values  the  criteria  (126)  and  (127)  are  easily  fulfilled  to  one  megahertz  and 
beyond,  making  very  attractive  the  consideration  of  a  compensator  for  film 
probes,  which  has  a  gain  proportional  to  Ju. 

11.  Conclusions 

The  main  conclusions  drawn  from  this  analysis  of  constant-current  film 
anemometers  are  as  follows: 

l)  In  contrast  to  hot-wire  probes.whlch  attenuate  the  input  as  the  simple 
Inverse  of  the  frequency  at  high-frequencies,  the  film  probe  attenuation  is  quite 
complex.  The  film  frequency  response  can  not  be  described  by  a  single  time 
constant;  rather,  this  response  depends  on  two  parameters:  the  inherent  film 
time  constant  and  a  loss  factor  dependent  on  the  substrate.  For  small  loss 
factors  (of  order  unity)  the  film  response  Is  almost  identical  to  the  wire 
response.  However,  the  loss  factor  for  actual  probes  is  estimated  to  be  of  order 
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100  or  higher.  In  such  cases  the  film  attenuates  more  than  the  wire  at  low 
frequencies,  although  its  response  falls  off  more  gradually,  at  a  3  db/octave 
rate  at  intermediate-to-high  frequencies.  At  extremely  high  frequencies, 
probably  high  enough  to  be  of  little  Interest,  all  films  resemble  the  wires  once 
again  as  regards  frequency  response. 

Two  practical  results  are  stressed.  First,  a  well -constructed  film  with 
small  inherent  time  constant  should  "see"  high  frequency  signals  quite  easily. 
Second,  since  the  ultra-high  frequencies  are  of  little  concern,  the  film  response 
of  3  db/octave  drop  may  be  easily  compensatable  by  an  electronic  circuit. 

2)  The  inherent  film  time  constant  resembles  functionally  that  of  the 
wire.  In  contrast  to  previous  analyses,  the  present  theory  shows  that  this  time 
constant  (for  the  wire  as  well)  should  include  a  term  related  to  the  dependence 
of  the  Nusselt  number  on  the  overheat. 

3)  The  phase  lag  angle  between  Input  and  output  is  identical  to  that  of 
hot  wires  for  a  loss  factor  of  zero,  and  constant  at  45  degrees  Independent  of 
frequency  when  the  loss  factor  Is  infinite. 

4)  The  sensitivity  coefficients  of  the  film  to  flow  fluctuations  (e.g.  in 
the  total  temperature,  velocity  etc.)  are  the  same  as  that  for  the  hot  wire. 
Forms  given  here  for  these  coefficients  alternative  to  those  in  the  classic  hot¬ 
wire  literature  are  presented  mainly  for  convenience. 

5)  The  steady-state  Nusselt  number  measured  with  the  film  by  the  classical 
overheat  method  Is  a  "compound”  one  which  Includes  conductive  losses  analogous 
to  the  "end  corrections"  for  hot  wires.  However,  whereas  these  wire  end 
corrections  are  numerically  inferior  to  the  actual  Nusselt  no.  measured  with  the 
wire,  in  the  case  of  the  film  they  can  be  substantially  higher  than  the  pure 
convective  Nusselt  number. 

12.  Recommendations  For  Further  Work 

In  order  to  embody  the  findings  of  this  report  into  a  practical  film 
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anemometer  system  capable  of  measuring  mean  and  turbulent  flows  at  high  speeds, 
the  following  extensions  of  the  work  done  to  date  are  recommended: 

1)  The  derivatives  ,  ffo,,  N„  etc.  appearing  in  the  operational  formulas 
(83),  (85)  etc.  must  be  evaluated  by  measurements  of  the  film  Mussel t  number  and 
recovery  factor  dependence  on  Mach  number,  Reynolds  number  apd  overheat.  The 
dependence  of  the  recovery  factor  on  Mach  number  is  judged  as  especially 
critical . 

2)  Once  these  quantities  are  known  over  a  range  of  the  flow  conditions, 
the  sensitivity  of  the  film  probe  to  flow  fluctuations  (e.g.  velocity  or  density 
fluctuations)  must  be  calculated  for  the  same  range  of  conditions. 

3)  The  calculated  3  db/octave  film  response  must  be  verified 
experimentally.  Compensating  amplifier  circuits  must  be  built  with  such  a 
3db/octave  gain,  to  replace  present  electronics  which  commonly  amplify  at  the  6 
db/octave  gain  characteristic  of  hot-wire  anemometers. 
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Temperature 


Figure  t.  Typical  hot-film  anemometer  probe  configuration  (above),  with  a  general 
view  of  the  temperature  distribution  in  the  probe  body  (below). 
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Plow  Conditions 
Temperature:  T0  +  T0’ (t) 

Nusselt  No.:  N  +  Nr  (t) 

Thermal  Conductivity:  ko  +  ko' 
Recovery  Temperature:  T*  =  qT0 

Flow 


Film  Conditions - 

Thickness:  Zero  (At  X  =  0) 
Temperature:  T  +  T'(t) 
Resistance:  R  =  Rr  [1  +  a  (T 
Heating  Current:  i 
Frontal  Area:  Af 
Heat  Capacity:  cf 


Lp 


■Back  Face 
(At  Constant 
Temperature, 
Te  =  iiT0) 


Substrate  Conditions 
Temperature:  Ts  (x)  +  Ts'  (x,t) 
Thermal  Conductivity:  ks 
Diffusivity :  (3 


Figure  2.  The  physical  model  with  definitions  and  nomenclature. 
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Hot-Film  Response 


Figure  3.  Attenuation  dependence  on  frequency  and  substrate  parameter  Q. 


Figure  4.  Phase  angle  dependence  on  frequency  and  substrate  parameter  Q. 
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Figure  5.  Example  of  a  film  probe  response  with  and  without  substrate,  and 
comparison  with  not-wire  probes. 
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LIST  OF  SYMBOLS 


Af:  Film  surface  area 

Ahu:  Hot-wire  surface  area 

A'„:  Morkavin's  overheat  parameter  (Reference  13) 

A:  Non-dimensional  amplitude  of  film  temperature  fluctuation 

B:  Non-dimensional  ohmic  heating  parameter  (eq.  (16}) 

c:  heat  capacity  (energy  per  degree  per  unit  mass) 

d:  Film  thickness 

dH:  Wire  diameter 

0:  Probe  diameter 

d(  ):  differential  (small)  fluctuation 

e:  Mean  film  voltage 

e':  Film  voltage  fluctuation 

e'„:  Net  film  voltage  fluctuation 

e, :  Film  sensitivity  coefficient 

E:  non-dimensional  input  amplitude 

f:  dimensionless  convection  transfer  parameter  (eq.  (15)) 

G:  Factor  expressing  sensitivity  to  fluctuations  (eq.  (78)) 

G:  Modified  form  of  G  (eq.  (80)) 

h:  film  height 

i:  electric  current 

I:  non-dimensional  current  fluctuation  amplitude  i'/i 

k:  thermal  conductivity 

1:  wire  length 

L:  penetration  length  of  film  heating  into  substrate  (eq.  (28)) 

Lp:  film  probe  length 

Lc:  A  characteristic  length 

m:  Exponent  In  temperature  dependence  of  conductivity 

M:  Mach  number 
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n:  Exponent  in  the  temperature  dependence  of  viscosity 

N:  Mussel t  number 

N:  compound  Nusselt  no.  {includes  conduction  losses) 

N„:  logarithmic  derivative  of  N  relative  to  M  (eq.  (61)) 

Np:  Nusselt  no.  for  convection  on  the  probe  body  surface 

Nr:  logarithmic  derivative  of  N  relative  to  r  (eq.  (61)) 

If,,,:  logarithmic  derivative  of  N  relative  to  Re  (eq.  (61)) 

q:  loss  factor  (eq.  (51),  (92)) 

r:  overheat  parameter  -  (R-R.)/Re 

?:  modified  overheat  parameter  (eq.  (81)) 

rp:  radius  of  probe  body 

r„:  material  electrical  resistivity 

R:  film  resistance 

Rr:  Reference  film  resistance  (at  1=0,  and  0°C) 

Re:  Reynolds  number,  based  on  stagnation  viscosity 

Re':  unit  Reynolds  number;  also,  fluctuation  in  Re 

Reh:  Reynolds  number  based  on  film  height  h 

t:  time 

t:  dimensionless  time  (eq.  (9)) 

T:  Temperature  (of  film,  if  unsubscribed) 

Tr:  film  reference  temperature  (taken  =  0) 

u:  flow  velocity 

w:  film  width 

x:  distance  from  film.  In  substrate 

a:  temperature  coefficient  of  resistance 

6:  substrate  material  dlffuslvlty  (-k/c p) 

& :  function  of  the  flow  Mach  no.  (eq.  (64)) 

e:  film  total  heat  capacity  (enetgy  per  degree) 

n:  temperature  recovery  factor 

tjm:  logarithmic  derivative  of  r\  relative  to  M  (eq. (61) ) 
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fjHt:  logarithmic  derivative  of  rj  relative  to  Re  (eq.  (61)) 

d  :  non-dimensional  temperature  (of  film,  if  unsubscribed)  see  eq.  (12) 

/i:  viscosity 

p:  density 

r:  Intrinsic  time  constant  of  film  (eq.  (41)) 

f:  non-dimensional  r  (eq.  (37),  (39)) 

phase  angle 

u:  angular  frequency 

u:  non-dimensional  angular  frequency  *  wG/Af 

(  )a:  at  recovery  temperature  conditions 
(  )(:  pertaining  to  the  film 
(  }D:  at  stagnation  conditions 
(  )4:  pertaining  to  the  substrate 
(  )w:  pertaining  to  a  hot  wire 

(  )':  fluctuating  quantity  (same  as  d(  ))  (note  Re'  exception). 

A 

(  }:  non-dimensional  fluctuation 
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